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SUMMARY 


An analytical investigation is made of the effect of shear lag on 
the bending vibration of vings that are designed essenti ally as shallow 
box beams, and a procedure is outlined for incorporating this effect 
in the determination of bending modes and frequencies. Numerical 
examples show that shear-lag action in a box beam can have a large 
influence on its vibration characteristics. The calculations indicate 
that even though only a small shear-lag action may be observed in a 
simple static deflection test of the beam, reductions in the second and 
higher-mode frequencies may be relatively large. 


INTRODUCTION 


Vibration tests of airplane wings have shown that discrepancies often 
exist between observed and calculated natural frequencies of wings. Among 
the possible sources of these discrepancies are aerodynamic and structural 
damping, rotary inertia, shearing deflections, and shear-lag effects, all 
of which are neglected in the usual engineering theory for beam vibration. 
The present paper investigates briefly the shear-lag effects. The shear- 
lag theory upon which this paper is based is presented in references 1 
and 2. 

The strength element of many wings is essentially a shallow box 
beam in which the secondary strains arising from shear lag sometimes have 
a significant influence on the bending stiffness, which in turn affects 
the vibration characteristics of the wing. That the natural frequencies 
of wings can be appreciably reduced when shear lag Is present is shown 
by the included numerical examples. A simple procedure is outlined for 
incorporating the effects of shear lag on bending stiffness in the determi- 
nation of bending modes and frequencies. Although not presented here, a 
similar procedure can be used to Incorporate the effect of the so-called 
bending stresses due to torsion on the torsional modes and frequencies 
of box beams. 
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K,p,q,7 


length of team. 

Young's modulus of elasticity 

effective shear modulus of cover sheet of toz team 
tending moment of inertia of team cross section 
mass of team per unit length 

equivalent loading used in numerical integration 
total shear force 
tending moment 

deflection of nth mode of vitration after i iterations for 
toth derived values of deflection and for values written in 
terms of unit tip deflection 

stations along team 

circular frequency of nth natural mode of vitration, radians 
per second 

frequency of nth natural mode of vibration, cycles per second 

number of equal-length bays into which team is divided for the 
shear -lag and vibration analysis 

distance between stations dividing toz team into bays 

statically indeterminate forces in corner flanges at each 
station due to shear lag 

radius of curvature of elastic axis at each station due to 
shear-lag strain in one cover of toz team 

parameters used in shear-lag analysis 


Ap area of corner angle or flange plus one-sixth the area of 

vertical shear wet 


A^ area of longitudinals and effective sheet material over one-hal 

the width of box (See reference !•) 
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b width from, center line of section to centroid of corner flange 

h depth "between centroids of top and bottom, flanges 

t average thickness of cover sheet 


MEANING OF SHEAR LAG 


Before an attempt is made to describe the effect of shear lag on 
the vibration of beams , a brief explanation of the term "shear lag" will 
be given. In figure 1 the tensile-stress distribution at two locations 
is shown for a simple sheet -stringer panel under a given loading. In 
the region near the load, the tensile stress is large at the edges where 
the load Is applied, but in regions remote from, the load the tensile 
stress tends to become more uniformly distributed over the width of the 
panel because of the shearing stresses which also exist in the structure. 

If these shearing stresses caused no shearing deformation of the Sheet 
in the panel, the tensile-stress distribution would necessarily be 
uniform at all points instead of varying as shown in figure 1. Shear lag 
is the term commonly used to describe the Influence that shearing deformations 
have on the stress distribution. 

The extent to which shear lag occurs in a structure Is a function of 
the geometry of the structure and of the loading. In the elementary 
bending theory of beams, the influence of shearing deformations on the 
stress distribution Is neglected because shearing deformations are generally 
small. In box beams of certain proportions, however, the shearing 
deformations cannot be Ignored if the stresses and deflection are to be 
predicted accurately.' A typical bending -stress variation over the cover 
sheet in the vicinity of the root of a cantilever box beam is shown in 
figure 2. Instead of a uniform longitudinal stress across the section, 
as would be predicted by elementary bending theory which considers the 
thin cover sheet infinitely stiff in shear, the longitudinal stresses are 
Increased near the edges and are decreased at the center of the beam. 

This change in stress distribution causes a change in the beam deflection. 

If cut-outs or certain load concentrations were present, similar shear- 
lag effects would occur. The deformation of wing structures, including 
the effect of shear lag, is discussed in reference 2 . 


EFFECTS OF SHEAR LAG ON VIBRATION 


Beam deflections are usually considered to be a function of the 
loading, the manner of support, and the beam stiffness which is calculated 
from the geometry of the cross sections and the modulus of elasticity of 
the material. The stiffness of box beams, however, is also Influenced by 
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shearing deformation of the thin cover sheets* For "beams in which this 
shearing deformation is appreciable, the increases in stress and, 
consequently, in strain in the corner flanges result in larger deflections 
for the beam as a whole than would be predicted by elementary beam theory. 

Such beams are essentially less stiff than similar beams (same El variation) 
without shear lag. In problems dealing with shear-lag beams, therefore, 
the usual concept of beam stiffness must be replaced by an effective 
stiffness concept, which takes into account the shear-lag strains present. 

With such a concept, the effective stiffness is a function of the beam 
loading, since the amount and distribution of the shear-lag strains along 
the beam vary for each different loading condition. 

Under the inertia loads occurring in vibration, the effective stiffness 
often differs appreciably from the geometrical stiffness. Since the stiffness 
characteristics of the beam change the vibration characteristics also 
change. Thus, shear-lag effects cannot be neglected in the determination 
of the bending modes and frequencies of same box beams. 


ANALYTICAL PROCEDURE FOR INCLUDING SHEAR -LAG EFFECTS 
IN THE DETERMINATION OF NATURAL MODES 
AND FREQUENCIES OF WINGS 


When bending modes and frequencies of beams are determined by the 
iteration procedure given in reference 3, the vertical shear forces and 
bending moments in the beam are found by direct Integration of an Inertia 

loading, m(x)Y^°^(x), where y( 0 )(x) Is the assumed deflection. The 
bending moments are converted into curvature by dividing by the stiffness Elj 
the slope and deflection are then found by integration of the curvature. 

This process is repeated, the newly-found deflection being used to compute 
the Inertia loading for the next iteration, .until two successively computed 
deflections bear a constant ratio to each other. This ratio is the square 
of the natural frequency o> . 

For a box beam with shear lag, the process is essentially the same 
except that the curvature cannot be found simply by dividing the bending- 
moment variation by the "geometrical" stiffness variation along the beam. 

M 

To the basic -“-curvature, a correction must be applied which takes into 
El 

account the additional curvature caused by the secondary strains in the 
corner flanges due to shear lag. The addition of this curvature correction 

to the basic ~-curvature is simply the process of properly taking into 
El 

account the effective reduction in stiffness of the beam. The resulting 
total curvature would be the same as if the bending moment had been divided 
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by the effective stiffness- Since the stiffness term El that appears 
in the differential equation of vibration for a beam is the stiffness 
that governs the curvature, the analysis of reference 3 still applies to 
the beams considered herein. In reference 3 the geometrical stiffness El 
computed from the geometry of the cross sections of the beam vas used, 
whereas in the present paper an effective stiffness which incorporates the 
effect of shear lag is used- 

An outline of the procedure used in this paper to determine the 
shear-lag curvature corrections will be given to clarify the numerical 
examples which are presented. The procedure is an adaptation of the more 
general shear-lag analysis presented in references 1 and 2. For a doubly 
symmetrical box beam, the steps are as follows: 

(1) Divide the box into 3 equal bays. (See fig. 3 .) For each 
bay, determine the constants Af, Al, b, h, and t for the simplified 
cross section (fig- 4-), as used in the shear-lag analysis presented in 
reference 1 . 

(2) From these simplified cross sections, compute the following 
shear -lag parameters: 

*-&Cr + £) 

Eb V A F k lJ 



t tanh EX 


q. 


K 

t sinh KX 


7 



(in reference 1, p, q, and 7 are defined with a factor ' G in the 
denominator; but, since G is a constant that can be factored from the 
equations given in step (4), the factor is omitted in these formulas.) 

( 3 ) From an assumed inertia loading mSf^ on the beam in the first 
iteration of the vibration analysis compute the average vertical shear 
force existing in the webs in each bay (one -half the total vertical 
shear for a given bay) . 
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(4) From the average vertical shear in the webs of each bay compute 
the self-equilibrated, groups of X-forces at each station. They are 
defined by the equations from reference 1 


X J 4 J " X J-1 (f.3 + + X ^-2lj-l 




“ + 7 J-1 


X J~lVl ‘ X J-2 (jj-l + P J-{) + X J-3<lj-2 = -7J-1 + 7 J-2 


>(D 


x l^l " X 0Pl 


= “ 7 ! 


J 


The subscripts on the X-forces refer to the stations and the subscripts 
on the parameters p, q, and 7 refer to the bays as they are numbered 
in figure 3- 

(5) Compute the curvature correction at each station from the 
X-forces by the relation 


1 _ x 
R “ EAj'h 


( 2 ) 


where Ajr and h are computed at the station points. Equation (2) 
gives the correction due to shear lag in one cover. 

(6) To the -curvature at each station computed in the vibration 

analysis, add the shear-lag curvature correction for both covers. Complete 
the remaining steps of an iteration. 

New curvature corrections should be computed in each succeeding 
iteration^ but if the inertia loading computed from the assumed deflection 
in the first iteration is reasonably representative of the true loading, 
the curvature corrections computed from this loading are accurate enough 
in most cases. Furthermore, precise computations are not 'justified because 
the accuracy of the shear-lag theory for loadings of the type that occur 
in vibration is unknown. In those cases for which curvature corrections 
are computed from assumed deflections differing widely from the derived 
shape, however, a second set of curvature corrections may improve the 
accuracy of the derived mode and frequency. 
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NUMERICAL EXAMPLES 


In order to show the effect of shear lag on "bending vibration, 
frequency and mode calculations have been performed for the two uniform, 
box beams shown in figures 5 and 6. Beam A (fig. 5 ) was designed to 
show only a smal l, amount of shear-lag action in static loading as shown 
by the fact that under a concentrated tip load the calculated tip 
deflection of this beam is only two percent greater by shear-lag theory 
than by ordinary beam theory. Beam B (fig. 6) was designed to show an 
exaggerated shear-lag action in static loading tests. Under a concentrated 
tip load, the calculated increase in tip deflection due to shear lag in 
this beam is about 19 percent. The following table gives a comparison 
of same of the natural bending frequencies of the two beams computed 
both with and without talcing shear lag into account: 


Beam 

Mode 

Frequency 

(cps) 

Percent 

change 

Shear lag 

No shear lag 


1 

1 ^ 5*7 

h 6.65 

2.1 

A 

2 

26b 

292 

12.3 


3 

672 

817 

21.6 

B 

1 

3 h. 6 

43.8 

21.0 


The calculations Indicate that, even though a box beam shows little 
shear-lag action under a simple static loading (for example, beam A), 
appreciable reductions in the second and higher-mode frequencies of vibration 
might be expected. Eor a box beam that shows large shear-lag action un der 
a simple static loading (for example, beamB), a correspondingly large 
reduction in even the fundamental bending frequency might be expected. A 
brief discussion of the calculation for each beam is presented in the 
following sections. 


Beam A 

Fundamental mode . - The procedure followed in the determination of 
the fundamental mode and frequency is indicated In the upper half of table 1 - 
The first iteration is shown in detail and each step follows closely the 
equivalent -load method shown in table 1 of reference 3, except for the 
insertion of the shear-lag curvature corrections (column 7), which are 
added to the -j-curvatures in column 8. (The mass per unit length m is 
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carried as a common factor in each column because it is constant over the 
length of the beam.) The shear-lag curvature corrections are calculated 
from the X-forces which are defined by the system of equations in the 
lower half of the table. The shear data in column 4 are used to compute 
the values of the parameter 7, which form the constant terms in the 
equations (equations (l))« The coefficients of the X-forces are determined 
from the shear-lag parameters p and q. The X-forces are solved for 
from this system of equations. The curvature corrections E/R are then 
calculated frcm these X-forces by the use of the equation shown in the 
table. (12/mX. 2 is a factor used to convert the curvature corrections to 

the same units as the -^-curvature and the factor 2 takes care of the 
correction for both covers of the symmetrical box beam. ) The curvature 
corrections are then inserted in the iteration procedure (column 7)* 

The fundamental mode, as determined from one more iteration, 1 b 
listed in column 15* Shear lag has little effect on the shape of the 
fundamental mode, as shown in figure 7(a) • The frequency of vibration 
with shear lag taken into account is computed in the table. If no shear 
lag were present, the fundamental frequency of this beam would be the 
same as that for a uniform cantilever, or 



3.52 


2 El 
mL 


4 


o 

= 85,700 (radians/sec) 


Second and third modes . - The procedure for finding the second and 
third modes and frequencies is the same as that given for higher-mode 
determination in reference 3 except for the shear-lag curvature corrections, 
which are introduced in the same manner as illustrated for the first 
mode in table 1. The computations show that shear lag accounted for about 39 
and 47 percent of the total curvature at the root station of the beam in 
the second and third modes of vibration, respectively. The curvature at 
the other stations was affected to a lesser but still significant extent- 
The effect of these curvature corrections on mode shape is illustrated 
In figure 7> It is evident that conclusions regarding the frequency 
change associated with these mode shapes cannot readily be drawn from 
comparison of the deflections alone because the curvature differences in 
the mode shapes, which also influence the frequency, are hidden. 
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Beam B . 

The length-width ratio, as well as the Af/A l ratio, for beam B 
is considerably smaller than, that for beam A and, as might be expected 
from the static-load deflection comparison presented previously for these 
two beams, a more pronounced shear-lag action was noticed for beam B 
in the fundamental -mode calculations. The calculations showed that shear- 
lag action accounted for about 48 percent of the total root curvature. 

The resulting sizeable reduction in frequency (21 percent) would not be 
expected in a box beam of ordinary design, but this reduction does 
indicate that the effect of shear lag can be appreciable even on the 
fundamental mode of vibration. 


CONCLUDING DISCUSSION 


The investigation has shown that, even though the shear lag in a 
given box beam had a relatively insignificant effect on the static 
deflection of the beam under a tip load and on the fundamental frequency 
of vibration, it caused an appreciable reduction in the higher frequencies. 
This increasing effect of shear lag on the higher modes can be explained 
by the fact that the relatively greater rates of change of the bending 
moment over the span in these modes of vibration are accompanied by 
increased shear deformation of the cover sheets and result in a reduced 
flexural stiffness of the beam. It has been pointed out that shearing 
deformations in the thin cover sheets effectively change the stiffness, 
and hence the vibration characteristics, of the beam- 

Shear deformation of the cover also takes place around discontinuities 
and abrupt changes in cross sections, such as cut-outs, and around points 
of load concentration. In an actual wing, shear-lag effects due to these 
disturbances may be of more importance than the simple effect treated in 
this paper and should therefore be investigated when a determination of 
the modes and frequencies is being made. The analytical procedure for 
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introducing shear-lag corrections for more complicated structures differs 
little from the procedure presented herein. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , February 2, 19^8 
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Figure 1.- Tensile -stress distribution in a panel with shear 

deformation of the sheet. 



Figure 2.- Bending-stress distribution in a cantilevered 

box beam. 
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Figure 3.- Notation for bays and stations in shear -lag and 

vibration analysis. 



Figure 4.- Simplified cross section used in shear-lag analysis. 
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Figure 7.- Effect of shear lag on the bending 

of a uniform box beam. 
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